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In supersymmetric seesaw model, the large flavor mixings of 
sleptons induce the lepton flavor violating (LFV) interactions 
iiljV (V — 7, Z), which give rise to various LFV processes. 
In this work we examine the LFV decays Z — > £i£j. Sub- 
ject to the constraints from the existing neutrino oscillation 
data and the experimental bounds on the decays £j — > lij, 
these LFV Z-decays are found to be sizable, among which 
the largest-rate channel Z — > tJX can occur with a branching 
ratio of 10~ 8 and may be accessible at the LHC or GigaZ 
experiment. 
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I. INTRODUCTION 

It is well known that the Standard Model (SM) pre- 
dicts an unobservably small branching ratio for any lep- 
ton flavor violating (LFV) process, such as £j — > £7-7 
or Z — > £i£j. In some extensions of the SM the LFV 
processes may be significantly enhanced [1-3]. One ex- 
ample of these extensions is the popular weak-scale su- 
persymmetry (SUSY). In SUSY models the LFV inter- 
actions IrljV (V = 7, Z) [4-7] receive two kinds of ad- 
ditional loop contributions: One is from the charged- 
current lepton-sneutrino-chargino couplings; the other 
is from the flavor mixings of charged sleptons. While 
the former is a common feature of all SUSY models ac- 
commodating right-handed neutrinos, the latter is siz- 
able only in some specific realizations of SUSY, such 
as the minimal supergravity model (mSUGRA) [8] with 
seesaw mechanism to generate the tiny masses for light 
neutrinos. The mechanism is realized by introducing 
right-handed neutrino superfields [1,2] with very heavy 
Majorana masses. In such a framework the flavor di- 
agonality of charged sleptons is usually assumed at the 
Planck scale, but the flavor mixings at weak scale are 
inevitably generated through renormalization equations 
since there is no symmetry to protect the flavor diagonal- 
ity Such flavor mixings of charged sleptons generated at 
weak scale are proportional to neutrino Yukawa coupling, 
which may be as large as top quark Yukawa coupling due 
to seesaw mechanism, and are enhanced by a large factor 
log(Mp/.M 2 ) (M P is Planck scale and M. is the neutrino 
Majorana mass). Therefore, the popular mSUGRA with 
seesaw mechanism predicts large flavor mixings of slep- 



tons at weak scale, which will reveal their effects through 
some LFV processes in collider experiments. 

The aim of this article is to examine the LFV Z- 
decays Z — > lilj induced by slepton flavor mixings in 
the mSUGRA seesaw model. Given the possibility of the 
extremely accurate measurement of Z-decays in future 
experiments, the decays Z — > £j£j may serve as a sensi- 
tive probe for such a new physics model. 

We will use the existing neutrino oscillation data and 
the experimental bounds on the decay lj — » £77 to 
constrain the model parameters, and then evaluate the 
branching ratios of Z — > £\£j. We find that, subject to 
the current constraints, the channel Z — > t/2 can occur 
with a branching ratio as large as 10 -8 and thus may be 
accessible at LHC [9] or the GigaZ option of TESLA at 
DESY [10]. 

This article is organized as follows. In section II, we 
briefly describe the SUSY seesaw model with minimal 
CP-violation in the right-hand neutrino sector and dis- 
cuss the induced flavor mixings between sleptons. In sec- 
tion III, we present the analytic results for the SUSY 
contributions to the branching ratio of Z — > liij. In 
section IV, we present the correlation between the pro- 
cess Z — » £j£j and lj — > £]j. In section V, we evaluate 
the numerical size of the branching ratio of Z — > lilj. 
Finally in section VI, we give our conclusion. 

II. SUPERSYMMETRIC SEESAW MODEL AND 
CHARGED SLEPTON MIXINGS 

A. Supersymmetric seesaw model 

The seesaw mechanism [11] provides an elegant expla- 
nation for the tiny masses of light neutrinos, which im- 
plies that new physics scale is about 10 14 GeV. However, 
a non-symmetric seesaw model suffers from a serious hi- 
erarchy problem [1], which can be automatically solved 
in the SUSY framework. 

In supersymmetric seesaw model with N right-handed 
neutrino singlet fields vr, additional terms in the super- 
potential arise [1]: 

W v = -\v£Wlv c R + Vr^ v L ■ H 2 , (1) 

where M is N x N mass matrix for the right-handed 
neutrino, and L and H2 denote the left-handed lepton 
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and the Higgs doublet with hypercharge —1 and +1, re- 
spectively. At energies much below the mass scale of the 
right-handed neutrinos, the superpotential leads to the 
following mass matrix for the left-handed neutrinos: 



M„ = mlM^mo = Y^M^Y^v sin/3) 2 



(2) 



Obviously, the neutrino masses tend to be light if the 
mass scale M. of the matrix M is much larger than 
the scale of the Dirac mass matrix = Yj,(iJ2) = 
Y„wsin/3 with v = 174 GeV and tan/3 = (H°)/{H°). 
The matrix M„ can be diagonalized by the MNS matrix 
U„: 

UtM^U* = diag(m vl ,m V 2, ra„ 3 ) , (3) 
where m V i are the light neutrino masses. 

B. Slepton flavor mixings 

The mass matrix of the charged sleptons is given by 

(4) 



m? m 2 ^ 
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m| + (to 2 - m 2 z sl v cos 2/3) 1, 
AfUCOs/3— m^/x tan /3 1, 



(5) 
(6) 
(7) 



where s^/ = sin 9w, Cw — cos 9w, 6w is the Weinberg 
angle and 1 is unit 3x3 matrix in generation space. In 
mSUGRA model it is assumed that at the Planck scale 
the soft-breaking parameters satisfy 



m, 



m„ — mz — to 1, m Hl = m H2 



mo, 



A e = A Y e , A V = A Y V . 



(8) 



In general, the lepton Yukawa couplings Yi and Y„ can- 
not be diagonalized simultaneously. It is usually assumed 
that Y^ is flavor diagonal but Y„ is not. In this basis 
the mass matrix of the charged sleptons is flavor diagonal 
at Planck scale. However, when evolving down through 
renormalization group (RG) equations (see Appendix A) 
to weak scale, such flavor diagonality is broken. In the 
leading-log approximation, we have [2] 



*K)" * -^(3m 2 + ^ 2 )(Y°tY°) / , 7 ln 
5{m%)u ~ , 
S(A e ) u : 



^(v?)„(Y?Y2)„m(^) 



(9) 
(10) 
(11) 



where Y° = Y(M P ). 



The flavor non-diagonal mass matrix m| in Eq.(4) at 
weak scale can be diagonalized by a unitary matrix S| 



(12) 



Such a unitary rotation of slepton fields is to induce 
the flavor-changing neutral-current vertices: Xa^i^x and 
Zlxly. 

In supersymmetric seesaw model, there exist right- 
handed sneutrinos with the same order masses as the 
heavy Majorana neutrinos. However, due to their large 
masses, they do not give significant contributions to 
the considered LFV processes. Therefore, only the left- 
handed sneutrinos need to be take into account, whose 
mass matrix is given by 



1 



m , 



to 2 ^ cos 2/3 1. 



(13) 



Due to the non-diagonal contribution (5(m|)/j in Eq.(9), 
m 2 is flavor non-diagonal at weak scale and needs to be 
diagonalized by a unitary matrix Sp 



S p m?St = diag(ml ). 



(14) 



Such a unitary rotation of sneutrino fields results in the 
charged-current flavor-changing vertex: X^tiVx- 



C. The form of neutrino Yukawa coupling 

As shown in Eqs.(9) and (11), the flavor mixings of 
charged sleptons are proportional to neutrino Yukawa 
couplings. Lack of knowledge of the neutrino Yukawa 
couplings results in numerous speculations on their pos- 
sible forms. Different forms may lead to different flavor 
mixings. In this work we consider a scenario called as 
minimal CP violating seesaw model which has two heavy 
Majorana neutrinos with the Dirac mass matrix m/j pa- 
rameterized as [12] 







ml = Yl{H%) = U L mV fi , m: 



TO 2 U 
TO 3 



where 
V fi = 



COS Or Sin Or 
— sin Or COS 8r 



e -i7R/2 o 

e^«/ 2 



(15) 



(16) 



with mixing angle Or and CP violating phase ^r for the 
heavy Majorana neutrinos concerning directly with lep- 
togenesis [12]. The matrix Ul appearing in Eq. (15) 
reads 

U L = 23 (6'L23)Ui2(6l I ,i3,5L)Oi2(e i i 2 )P(-7 I ,/2) , (17) 

where P(-7 L /2) = diag[l, cxp(-z7 L /2), exp(z7 L /2)], 
and Oij and denote the rotations in (i,j) plane. 
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Without loss of generality, m2,3 in Eq.(15) are chosen 
to be real, positive and m 2 < W3. 

As Eqs. (2) and (15) are used, the mass matrix for the 
light neutrinos in this model can be further expressed as 



M„ = UimVflM-'V^m 7 ^ 



(18) 



And the MNS matrix in (3) is found to be a product of 
matrices, 



(19) 



where Kj? : = 'K. R (9 1 <f>, a) is a unitary matrix. Therefore, 
now Eq. (3) can be rewritten as 



K^mV fl M- 1 V^m T K| ? = diag[ 



(20) 



From this equation, one can learn both K/{ and m Ui are 
independent of the choice of XJl- 

It is noticeable that the special form (15) for the neu- 
trino Yukawa couplings matrix ~Y V implies [12]: 

(1) One of the neutrinos is massless, i.e., m v \ = 0. 

(2) The quantity Y^Yj, is only dependent on 3 mixing 
angles 0li2,li3,L23 and a CP violating phase Sl in Ul, 



(YiY„ 



m§(U £ ) f2 (U* ) 2J + ml(y L ) I3 (Vlh 
(v sin /3) 2 



(21) 



(3) For small mixing angles 0li3 and 9, the light neutrinos 
mixing matrix U„ takes a simplified form similar to the 
mixing matrix introduced in [13]: 



U„ 



CL12 SL12 

— SL12CL23 CL12CL23 SL23 
\ S_L12S_L23 — CL12SL23 C£23 



A 

P(a') (22) 

2CL23 S£23 
12SL23 C_L23 / 

where <f>' — <fi + jl, a' = a — 7l/2 and s x = sinx, c x = 
cos a;. In this case, the angles in Y^Y^ can be related 
directly to the corresponding neutrino mixing angles and 
determined by neutrino experiments. 

III. THE LFV DECAYS Z -> liij 

The flavor changing interactions in slepton sector dis- 
cussed in the preceding section, namely the couplings 
Xa^i^J and Zljlj from charged slepton mixings as well 
as xi^i^j from sneutrino mixings, can induce the LFV 
processes Z — > £j£j, as shown in Fig. 1. The relevant 
Feynman rules can be derived straightforwardly from the 
analysis in the preceding section. Our analytic results 
will be expressed in terms of the constants g^g 1 , Gxy 



and Cr^t defined in Fig. 2, whose explicit expressions 
can be found in [7,14]. 




FIG. 1. Feynman diagrams of SUSY contributions to the 
LFV processes Z — > Zilj. 
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FIG. 2. Some interaction vertices needed to calculate the 
branching ratio of Z — > ijlj in SUSY, a and /3 are indices of 
charginos (neutralinos) , while X and Y are those for sleptons. 



The calculation of the diagrams in Fig. 1 results in an 
effective Zljlj vertex: 

M = ige^unj {p x )Yi M u ij (p 2 ) (23) 

with £ M being the polarization vector of Z-boson, p\{p2) 
the momentum of £i(ij), and r M given by 
a, 



sin 6w 



[7 M (AlPl + IirPr) 



ll/ k u (f 2L P L + f 2R P R )} , 



(24) 



where Pr,l = |(1±75), g = e/ s'mBw and k — pi — p 2 is 
the momentum transfer. The form factors fiL, fi R , f 2 L 
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and fin arising from the calculation of the loop diagrams 
in Fig. 1 are listed as follows. 
Contribution of Fig. 1(a): 



1L — ^XY^IaX L /0 24°JaY 

+m ej m a {C a v+C a 11+ C a 12 )C L }aY ], 
f%L = GxyC*j^ x [m a (Cq + Cl x + Ci 2 ) Cj aY 
-m e j (Ci 2 + C'22 + C 23 ) Cf aY ] ■ 

Contribution of Fig. 1(b): 
fL = CtxC% x [g^m a m C b 



(25) 
(26) 



m Z°23 



2C b 4 + - 



+C^ x C% x g^m a m ej (C b + C\ x + C h 12 ) 
fL = C** X C% X {g%m a C b 1+ g L ap m C b 2 ) 

+ ClaxCfp X 9a/3 m ej (Cl2 + C22 + C23) ' 

Contribution of Fig. 1(c) plus Fig. 1(d): 



1L — WaX 



mj 



-{B - tS Q )C JaX - tfiC JaX 



9L, 



. 



(27) 
(28) 

(29) 
(30) 



In the above, gh = (1 — 2 sin 2 #w)/(2cos#w), and Bq x = 
B(- Pl :m 2 a ,m x ), Cg^ = C ,ij (-pi, -pi\ m 2 , m Y , m\) 
and = Co ti j (-pi, -p 2 ; m 2 Xl m 2 f3 , to 2 ) are the Feyn- 

man loop integral functions [15]. Terms proportional to 
the lepton masses mt I are neglected. The right-handed 
form factors from the vertex loops are obtained from the 
corresponding left-handed ones in (25)- (28) by the sub- 
stitution L <-> R. 

The branching ratio of Z — > £j£j (including its charge- 
conjugate channel) is then given by 1 



Br(Z - = 



a e 



48tt 2 Vsin 2 ^/ T z 



m z 



lif- 



1L 



+^(\f2L\ 2 + \f2n\ 2 ) 



l/l«| 2 

(31) 



where fiL,%R = f?L iR an d denotes the total 

a— a,6,c 

decay width of Z boson. 

Although the above results are sufficient to allow for 
numerical calculations, we would like to derive an ana- 
lytical expression for the branching ratio by considering 
the limit ms mz where ms represents the mass of 
any internal sparticle in the loops in Fig.l. In this case 



1 Our result is in agreement with that given in [7] if mtj 
dependence terms in fiL,iR are neglected. 



the loop functions can be much simplified and we use the 
mass-insertion approximation in our derivation. In such 
a limit, the chargino mass matrix 



M2 \f2mw sin (3 
\/2m\v cos (3 yu 



(32) 



is nearly diagonal. Here fi is the mass parameter appear- 
ing in the term [j,HiH 2 in superpotential and M 2 is the 
SU(2) gaugino mass parameter. The matrices U and V 
which diagonalize M^± will be unit ones for \x > 0, and 
the chargino masses are given by 



m-± = M 2 , m-± 

Al X2 



(33) 



The symmetric neutralino mass matrix 
/ 



Mi 

M 2 

mzswc/3 mzcwc/3 

y mzswsp -mzCwS/3 -// J 



(34) 



can be diagonalized by a unitary matrix N 



N = 



/I 



V 



V2e^ 
-V2V 4 ' 



V2e l f J 



The corresponding mass eigenvalues are given by 

Ml i2 , m-fO = TTLyO = \p\. 



x.1,2 



(35) 



(36) 



When using mass-insertion method, one should note 
the fact that, for any matrix M = M° + M 1 , where 



M° = diag(TO?, 



and M 1 has no diagonal el- 



ements, if matrix T can diagonalize the matrix M, 
TMT 1 = diag(mi,TO2, ■ ■ ■ ,to„), then at leading order 
for an arbitrary function / 



T\ k f(m k )T k0 = 5 tf /(m?) + M^(m°,m°) 
with 

j{x,y,Zf-z n ) = 

x-y 



(37) 



(38) 



After a straightforward calculation we obtain an analyt- 
ical expression for the branching ratio 



Br(Z -» £dj) 



aim c 2 w m z \S(ml)i 
M| 



48tt 2 s 6 w Tz 



\fi(xi,xj) 



-2f 2 (xi,xj) 



12 A 
2 S W — S W 



5 + Sw h( x .y XT ) 



M 

Mi 



CjX'j 



„2 



h{xi,xj) 



(39) 
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Here Sw — sin9w, c\v — cos#w, xj — (m~)///M|, 
x'j = (m|)///Mf, and 



h{x) = —^—r { 1 ^-r In r 



x — 1 V x— 1 



1 / x 2 - 2x , 
/ 3 (a:) - - I 1 + — In a; 



(x- 1) V a; - 1 
and fi(x,y) can be obtained through Eq. (38). 



(40) 



IV. COMPARISON OF LFV Z-DECAYS WITH 
LEPTON DECAYS 

Now we compare the LFV Z-decays with lepton de- 
cays. Using a similar procedure in the preceding section, 
we can easily calculate the decay width for tj — > £jj by 
setting g = e, g^'J 1 = 1 with a = (3, Gxy = 1 with 
X = Y in Fig. 2, and f 1L . 1R = in (24). Meanwhile 
one should also note the fact that sneutrinos in Fig. 1(a) 
and neutralinos in 1 (b) do not couple to photon and that 
the self-energy diagrams do not contribute to dipole op- 
erators. The branching ratios of £j £i-f are obtained 
as 



Br(£j -> 1/7) _ 6a em .. , 
Br(£j — > ijvjvi) 7r m 2 



(l/ 2 7 J 2 + l/ 2 V)- (41) 



Here the form factors are given by [7] 

- E ^ C "x [r^C^Ft + m tj C«&F*\ , (42) 

1 L L J 

fe — a,t> 

= / 2 \ \l»r , (43) 
where 



F 1 a (x a )=m^o(C a + C 1 a 1 +C7 1 a 2 ) 



lna; a ) , 



(l. — l) 2 V 2 Xa - 1 

F 2 a (x a ) = -m?o (Cf 2 + C 2 a 2 + C 2 a 3 ) 

~X a ~t~ 5x'q, ~\~ 2 



(44) 



2(:r a - l) 3 

. 2 



F?(x b )=m^(C b 11 + C b 12 ) 



— 3a; b + 1 x 



(x 6 - l) 2 V 2 x b - 1 



£a - 1 



In x b , 



lnx a , (45) 



(46) 



F\ (x b ) = m 2 (C b 12 + C b 22 + C b 23 ) = - F 2 a ( - ) (47) 

Xcv Xb Xb 

with x a = m?^/m| and x 6 = m? x /m|_. 

Next we derive the analytical expression for the 
branching ratios in the limit of mj >> mz- Unlike the 
form factors for the Z-decays which contain terms not 
proportional to the small lepton mass [see Eqs. (25) and 



(28)], the form factors for lj — ► £jj are always propor- 
tional to the small lepton mass mi,. In this case, the 
off-diagonal elements in the mass matrices of chargino 
and neutralino are no longer negligible, especially when 
tan/3 is large. In fact, the terms m a Cj^^ Cfax m fin, 
receive the contribution from the wino-Higgsino mixing, 
which can be enhanced by tan (3. So for a large tan /3, the 
contribution of f^ R is dominant and the branching ratios 
are given by [2] 

Br(£j -> £ n ) ~ Br(£j -» livjvi) 
= Br(£j -» l lVj vi) 



6a e m 


m w 


l/ 2 7 «| 2 


ir 


ml 




6a em 


m w 




ir 







-F?(x It xj)-Fi(xi,xj) 



\S{jn\)u\ 2 tan 2 /? 



F®(xi, xj) — F^(xi, xj) 



(48) 



(m|WM 2 - 



where xj 

Comparing Br(£j — ► £77) with Z5r(Z — > £i£j), we find: 

(1) The dipole transitions in (24), the only operators 
contributing to £j — ► £77, do not give dominate contri- 
butions to decays — > £/lj due to heavy sparticle mass 
suppression; 

(2) £>r(Z — > £/^j) is not sensitive to tan/3, whereas 
Br(£j — > £77) can be enhanced by large tan/3; 

(3) The ratio £r(Z -» tih)/Br{£j -> £77) is indepen- 
dent of the heavy Majorana sector introduced by seesaw 
mechanism. 



V. NUMERICAL RESULTS 

In our numerical calculation we consider the con- 
straints from current neutrino oscillation experiments 
and the experimental bounds on LFV lepton decays. 

(1) Neutrino oscillation experiments: 
The SK Collaboration [16] showed that the created 
in the atmosphere oscillates into v T with almost maxi- 
mal mixing, sin(2# atm ) <~ 1 and the neutrino mass-square 
difference is Aml tm ~ (2 - 4) x 1(T 3 eV 2 . The second 
mass-square difference and mixing angle are found to be 
Am L = ( 3 - 15 ) x 10_5 cV2 ; sin(26» so; ) = 0.7 ~ 0.9 
from solar neutrino experiments [17,18]. For the third 
mixing angle, only the upper bound is obtained from the 

reactor neutrino experiments [19,20]: sin 2 29 rea < 0.1 for 
2 



Am 2 tm ~ 3 x 10- 3 eV 



Although there exists a possibility that neutrino 
masses are quasi-degenerate, in this work we take the 
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normal mass order m v \ < m v2 < m„ 3 with values: 2 



m vl = 0, m v2 = JAm 2 sol , m v 



AmL. (49) 



The mixing angles are fixed to be 

^12 = ^=30°, 6 L23 = e atm = 45°. (50) 

Further, we restrict 6 L13 < 10°. Then (YjY^/j in Eq. 
(21) are given by 



(YtY„) 12 
(YtY„) 13 



X — m2 + an2^i 3 m3 , (51) 



4u 2 sin 2 /? ^ 2 
V2 / V3 



4i> 2 sin 2 (3 



m\ + sin 20 L i 3 m 2 , (52) 



23 



4i> 2 sin 2 /3 



(2m^ - m§) . 



(53) 



The dependence of the parameter Y£Y„ on CP phase Sl 
is very weak and thus has been neglected. 




500 rooo rsoo o 500 woo woo 2000 

m (GeV) m (GeV) 

FIG. 3. Branching ratios oi Z —* iitj and £j — > £77 versus 
the common scalar mass m . Other parameters are fixed to 
be mi/2 = 800 GeV, A = 0, tan/3 = 10, m 2 = 10 GeV, 
ms/m,2 — 30 and 9lva = 0. The dash line in (b) is the 
experimental upper bound on /i — ► e7- 

Experimental bounds on LFV lepton decays: 
LFV lepton decays have been searched in several experi- 
ments and the current bounds are given by [22-25] 



Br{/j,^ej) < 1.2 x 10" 



(54) 



2 In general, the impact of RG evolution on neutrino masses 
and mixing angles can be large, however, it is small for a 
hierarchy of neutrinos we chosen [21]. 



Br(r -> (e, < (2.7, 1.1) x 10~ 6 , (55) 

Br{Z r/2) < 1.2 x 10" 5 , (56) 
Sr(Z -> (/i, r)e) < (1.7, 9.8) x 10~ 6 . (57) 

In addition, the explanation of the observed lepton 
number asymmetry by seesaw mechanism gives a lower 
bound for heavy Majorana neutrinos A4\ > 10 11 GeV 
[12]. Taking into the constraints mentioned above 
and fixing the right-handed neutrino masses as M.\ = 
10 13 GeV, M 2 ^ 10 15 GeV, we solve the full RG equa- 
tions listed in Appendix A numerically based on the work 
of [26] , where the experimental bounds from b — > s-f and 

— 2 have been already taken into account. Although 
the processes Z — > iilj are closely correlate to £j — > £77 
and there is a quite stringent bound on fi — > e7, our 
numerical results show that there exists a scenario with 
7712 <SC m3 and a very small #1,13, in which a large branch- 
ing ratio for Z — > r/I is obtained. 

In Fig. 3 we show the branching ratios of Z — > ^/^j 
and £j — > £77 versus the common scalar mass mo- From 
Fig. 3 we have the following observations: 

(1) With fixed m\i 2 and tan/3, both Br(Z — > r/i) and 
Br(r — » ^7) reach their maximum values as mo ~ 
1000 GeV, then drop slowly as m gets larger. 



(2) The branching ratio of Z 
HT 8 . 



Tfj, can be as large as 



Since 5.5 x 10 9 Z-bosons will be produced at the LHC 
[9] and the possible sensitivity of GigaZ to Z — ► Tp will 
be up to 10~ 8 [10], the mode Z — ^ Tp will be accessible 
at both the LHC and TESLA GigaZ. 




2.5 5 7.5 2.5 5 7.5 10 

9 LJ3 (degree) 9 LJ3 (degree) 

FIG. 4. Same as Fig. 3, but versus the mixing angle Olva 
with m = 500 GeV. 

It is noticeable that the branching ratios are sensitive 
to the mixing angle 9 lis except for the processes Z — ► Tp 
and r — > /z7. As an illustration, we plot the dependence 
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on 9li3 in Fig. 4. We see that to satisfy the experimental 
constraint on /i — ► ej, the mixing angle 9li3 must be 
quite small. Therefore, a joint measurements for LFV Z- 
decays and lepton decays will set strong constraints on 
the model parameter space. 



VI. CONCLUSIONS 

We evaluated the lepton flavor violation Z decays in 
the framework of supersymmctric seesaw model at first 
time. Although different forms of neutrino couplings may 
lead to different size of LFV Z decays, we emphasize that 
it is important to study how large the rate for the LFV 
can be for some typical cases and analyze the possibility 
to observe Z — > £[£j in future experiments. From our 
calculation results we conclude that, subject to the con- 
straints from the existing neutrino oscillation data and 
the experimental bounds on the decays £j — ► £ij, the 
LFV Z-decays Z — > £j£j can still be sizable in super- 
symmetric seesaw model, among which the largest-rate 
channel Z — > t£l can occur with a branching ratio of 
10~ 8 and thus may be accessible at the LHC and GigaZ 
experiment. 
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APPENDIX A: RENORMALIZATION GROUP 
EQUATIONS IN SUSY SEESAW MODEL 



(2) For soft parameters: 
1 



dm| 
dt ~ 16tt2 

+2m 2 H2 Yt Y v + 2A* A„] , 
dm? _ 1 
~dT ~ 8^ 

+2m^ 2 Y, y Y,t + 2A„At] , 

1 



rx& Yt Y v + Yj Y v m 2 ~ + 2Yt m ? Y„ 



m%Y v Yt + Y„Yt m? + 2Y v m~Y\, 



dm 2 H2 



dt 



8tt 2 



Tr 



Yt (m| + m?+m^ 2 ) Y„ 

+ At A„] , 

= -L(2Y / YtA v + A / Yt Yv ), 



dt 



16tt 2 



(A4) 

(A5) 

(A6) 
(A7) 



— = ^ { [T 2 - 5 ? - 3 ff2 2 + 4Y„Yt] A„ + A„Y] Y, 



+ [2Tr(3Yt r A t/ + Y+ A„) + 5A„Yt 
-2 {glMt + 3glM 2 ) Y v + 2Y,Y]A,} 

(3) For neutrino masses [27]: 



Y„ 



(A8) 



^ = _L[M(Y,Ytf +Y.Y+M], 



(A9) 



dM v 



6tt 2 I 



dt l^ll 2T2 + ^ +Y ^ T \ 

+M V (Yt Y, + Y]Y £ - 2 5 2 - 6 ff2 2 ) } . (A10) 

Note that the above RG equations are valid for the run- 
ning from Mp to M. Below the scale M, the RG equa- 
tions are the same except that the couplings of the right- 
handed neutrinos do not appear. 



In this appendix we present additional contributions to 
the RG equations of some parameters in supersymmetric 
seesaw model due to non-zero neutrino interactions. The 
detailed description of these equations can be found in 
[1,2]. At one- loop level, the RG equations are given as 
follows. 

(1) For Yukawa couplings: 

(T 2 - g\ - 3g 2 2 + 3Yt Y v + yJy*) , (Al) 



dY v 


Y„ 


dt 


16tt 2 


dY e 


Y e 


dt 


16tt 2 


dYu 


Yu 


dt 


16tt 2 



Yt Y 

Tr(YtY„), 



(A2) 
(A3) 



where t = In \x r with fi r being the renormalization scale, 
and T 2 = Tr(3Y^Yc; + Yj Y v ). Y v is the Yukawa cou- 
pling matrix for up-type quarks, and g±, g 2 and 173 are 
the U(l)y, SU(2) and SU(3) gauge coupling constants, 
respectively. 
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